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Abstract 

We consider the cubic nonlinear Schrodinger equation in two space dimensions with an 
attractive potential. We study the asymptotic stability of the nonlinear bound states, 
i.e. periodic in time localized in space solutions. Our result shows that all solutions with 
small, localized in space initial data, converge to the set of bound states. Therefore, the 
center manifold in this problem is a global attractor. The proof hinges on dispersive 
estimates that we obtain for the non-autonomous, non-Hamiltonian, linearized dynamics 
around the bound states. 

1 Introduction 

In this paper we study the long time behavior of solutions of the cubic nonlinear Schrodinger 
equation (NLS) with potential in two space dimensions (2-d): 

idtu{t,x) = [-A.^ + V{x)]u + -f\ufu, t > 0, x G (1) 
u{0,x) = uo{x) (2) 

where 7 G M — {0}. The equation has important applications in statistical physics, optics 
and water waves. It describes certain limiting behavior of Bose-Einstein condensates |S1 E] 
and propagation of time harmonic waves in wave guides • In the latter, t plays the 

role of the coordinate along the axis of symmetry of the wave guide. 

It is well known that this nonlinear equation admits periodic in time, localized in space 
solutions (bound states or solitary waves). They can be obtained via both variational tech- 
niques m 123 1^ and bifurcation methods f^l 1^ , see also next section. Moreover the set 
of periodic solutions can be organized as a manifold (center manifold). Orbital stability of 
solitary waves, i.e. stability modulo the group of symmetries u 1— > e"*^M, was first proved in 
iniES], see also jHllinilSll- 

In this paper we are going to show that the center manifold is in fact a global attractor 
for all small, localized in space initial data. This means that the solution decomposes into a 
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modulation of periodic solutions (motion on the center manifold) and a part that decays in 
time via a dispersion mechanism (radiative part). For a precise statement of hypotheses and 
the result see Section |3] 

Asymptotic stability studies of solitary waves were initiated in the work of A. Soffer and 
M. I. Weinstein j2Sll2ni, see also [21 El IH 13 IH] • Center manifold analysis was introduced in 
[in], see also [22] • The techniques developed in these papers do not apply to our problem. 
Indeed the weaker — > L°° dispersion estimates for Schrodinger operators in 2-d, see fl24|) . 
compared to 3-d and higher, respectively lack of end point Strichartz estimates in d = 2, 
prevent the bootstrapping argument in [3 IHl EHl EE] , respectively [TT], from closing. The 
technique of virial theorem, used in [3 El E] to compensate for the weak dispersion in 1-d, 
would require at least a quintic nonlinearity in our 2-d case. Finally, in [221, the nonlinearity 
is localized in space, a feature not present in our case, which allows the author to completely 
avoid any L°° estimates. 

To overcome this difficulties we used Strichartz estimates, fixed point and interpolation 
techniques to carefully analyze the full, time dependent, non-Hamiltonian, linearized dynamics 
around solitary waves. We obtained dispersive estimates that are similar with the ones for 
the time independent, Hamiltonian Schrodinger operator, see section |3] Related results have 
been proved for the 1-d and 3-d case in (201 [13 [22] but their argument does not extend to the 
2-d case. We relied on these estimates to understand the nonlinear dynamics via perturbation 
techniques. We think that our estimates are also useful in approaching the dynamics around 
large 2-d solitary waves while the techniques that we develop may be used in lowering the 
power of nonlinearity needed for the asymptotic stability results in 1 and 3-d mentioned in 
the previous paragraph. 

Note that, in 3-d, the case of a center manifold formed by two distinct branches (ground 
state and excited state) has been analyzed. Under the assumption that the excited branch 
is sufficiently far away from the ground state one, in a series of papers [23 [23 [23 
the authors show asymptotic stability of the ground states with the exception of a finite 
dimensional manifold where the solution converges to excited states. We cannot extend 
such a result to our 2-d problem as of now. The reason is the slow convergence in time 
towards the center manifold, t~^~^ in 2-d compared to in 3-d. This prevents us from even 
analyzing the projected dynamics on a single branch center manifold, i.e. the evolution of 
one complex parameter describing the projection of the solution on the center manifold, and 
obtain, for example, convergence to a periodic orbit as in |3 [21] . However the evolution 
of this parameter, respectively two parameters in the presence of the excited branch, is given 
by an ordinary differential equation (ODE), respectively a system of two ODE's, and the 
contribution of most of the terms can be determined from our estimates, see the discussion in 
Section [31 We think it is only a matter of time until the remaining ones will be understood. 

The paper is organized as follows. In the next section we discuss previous results regarding 
the manifold of periodic solutions that we subsequently need. In section [21 we formulate and 
prove our main result. As we mentioned before the proof relies on certain estimates for 
the linear dynamics which we prove in section [^ We conclude with possible extensions and 
comments in section [S] 

Notations: H = -A + V; 

LP = {/ : M2 ^ C I / measurable and J^, \f{x)\Pdx < oo}, = {J^, |/(x)|*'dx)'/^ 
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denotes the standard norm in these spaces; 

< X >= (1 + and for cr e M, denotes the space with weight < x >^°', i.e. 

the space of functions f{x) such that < x f{x) are square integrable endowed with the 
norm ||/(x)||i2 = II < X /(x)||2; 

{f^9) = /k2 f{.x)g{x)dx is the scalar product in where / = the complex conjugate of 
the complex number /; 

Pf. is the projection on the continuous spectrum of H in L^; 

H" denote the Sobolev spaces of measurable functions having all distributional partial 
derivatives up to order in L^, || ■ ||//n denotes the standard norm in this spaces. 

2 Preliminaries. The center manifold. 

The center manifold is formed by the collection of periodic solutions for (^: 

UE{t,x) = e-''''Mx) (3) 
where E and ^ tpE G i^^(M^) satisfy the time independent equation: 

[-A + V]tPE + i\^e\^^e = E^Pe (4) 

Clearly the function constantly equal to zero is a solution of (jH) but (iii) in the following 
hypotheses on the potential V allows for a bifurcation with a nontrivial, one parameter family 
of solutions: 

(HI) Assume that 

(i) There exists C > and p > 3 such that: 

|V^(a;)| < C < X forallxeM^; 

(ii) is a regular point^ of the spectrum of the linear operator H = —A + V acting on L^; 

(iii) H acting on has exactly one negative eigenvalue Eq < with corresponding nor- 
malized eigenvector ipQ. It is well known that tpo{x) can be chosen strictly positive and 
exponentially decaying as |x| oo. 

Conditions (i)-(ii) guarantee the applicability of dispersive estimates of Murata fH] and Schlag 
[221 the Schrodinger group e~*^*, see section |3J In particular (i) implies the local well 
posedness in of the initial value problem (012)), see section |21 

Condition (iii) guarantees bifurcation of nontrivial solutions of (jH) from [Eq, iPq). In Section 
13 we discuss the possible effects of relaxing (iii) to allow for finitely many negative eigenvalues. 
We construct the center manifold by applying the standard bifurcation argument in Banach 
spaces jlH] for ^ at E = Eq. We follow [1^] and decompose the solution of (H)) in its projection 
onto the discrete and continuous part of the spectrum of H : 

ipE = ai^o + h, a= {ipo, ^e), h = Pc^e- 
^see [23 Definition 7] or = {0} in relation (3.1) in [El 
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Using the notations 

fp{a, h) = {ipo, \aipo + hl'^^aipo + h)), (5) 
fc{a, h) = Pclatpo + /ip(a^o + h), (6) 

and projecting (jH) onto ipo and its orthogonal complement = Range Pc we get: 

h = -j{H-Er'aa,h) (7) 
Eo-E = --fa~'fp{a,h). (8) 

Although we are using milder hypothesis on V the argument in the Appendix of ^Hj can be 
easily adapted to show that: 

T{E, a,h) = h + -f{H - Ey^f,{a, h) 

is a function from (-cx), 0) x C x Llr\H^ to Ll^H^ and J^(Eo, 0, 0) = 0, DhJ^{Eo, 0, 0) = /. 
Therefore the implicit function theorem applies to equation ((7j) and leads to the existence of 
6i> and the function h{E, a) from {Eq - 6i, Eq + 6i) x {a e C : \a\ < 6i} to Ll n H"^ 
such that ((Tj) has a unique solution h = h{E, a) for all E G {Eq — Si, Eo + 6i) and \a\ < 6i. Note 
that if (a, h) solves © then (e'^a, e'^/i), 6 G [0, 27r) is also a solution, hence by uniqueness we 
have: 

h(E,a) = -^h(EAa\). (9) 

\a\ 

Because ipo is real valued, we could apply the implicit function theorem to ((Zj) under the 
restriction a G M and h in the subspace of real valued functions as it is actually done in jl9j . 
By uniqueness of the solution we deduce that h{E, \a\) is a real valued function. 

Replacing now h = h{E, a) in (jSJ and using (0) and Q we get the equivalent formulation: 

Eo- E = -^\ar'fp{\alh{E,\a\)). (10) 

To this we can apply again the implicit function theorem by observing that G{E, a) = Eq — 
E + ■y\a\~^ fp{\a\, h{E, \a\)) is a function jTHl Appendix] from {Eq — 5i, Eq + 6i) x (— 5i, 6i) 
to M with the properties G{Eo,0) = 0, dEG{Eo,0) = —1. We obtain the existence of < 
S < Si, < 6e < Si and the function E : (—6,6) i— > {Eq — 5e,Eq + 5e) such that, for 
\E — Eq\ < 6e, \a\ < 6, the unique solution of (jHI) with h = h{E, a), is given hj E = E{\a\). 
If we now define: 

h{a) = —h{E{\a\, \a\) 
\a\ 

we have the following center manifold result: 

Proposition 2.1 There exist Se,S > and the function 

h:{aeC : \a\ < 6} ^ LlnH^ 

such that for \E — Eq\ < 5e the eigenvalue problem ^ has a unique solution up to multipli- 
cation with e*^, 9 G [0,27r), which can he represented as: 

ipE = cii^o + h{a), {ipo, h{a)) = 0, \a\ < S. 
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Since ipo{x) is exponentially decaying as |x| oo the proposition implies that ipE ^ L"^- 
A regularity argument, see [23], gives a stronger result: 

Corollary 2.1 For any cr G M, there exists a finite constant C„ such that: 



We are now ready to prove our main result. 

3 Main Result. The collapse on the center manifold. 

Theorem 3.1 Assume that hypothesis (HI) is valid and fix a > 2. Then there exists an 
Eq > such that for all initial conditions uo{x) satisfying 



with the constants Ci,C 2, independent of Eq and C2 depending on p> 2. 

Before proving the theorem let us note that (fTT|) decomposes the evolution of the solution 
of (CQ)-© into an evolution on a center manifold ipE{t) and the "distance" from the center 
manifold r{t). The estimates on the latter show collapse of solution onto the center manifold. 
The evolution on the center manifold is determined by equation (fT^ below. We discuss it in 
Sectional 

Proof of Theorem 13. IL It is well known that under hypothesis (Hl)(i) the initial value 
problem (PJ-Q is locally well posed in the energy space and its norm is conserved, see 
for example Corollary 4.3.3. at p. 92]. Global well posedness follows via energy estimates 
from ||uo||2 small, see Remark 6.1.3 at p. 165]. 

In particular we can define 



max{||no||L2, ||no||//i} < £0 



the initial value problem (Op-IB' ^■^ globally well-posed in H^. 
Moreover, for all t and p > 2, we have that: 

u{t, x) = a{t)ipo{x) + h{a{t)) +r(t, x) 



(11) 




a{t) = {tpo,u{t)), hiteR. 



Cauchy-Schwarz inequality implies 



'oil 2 — 



\\uoh < ^0, for t G M 
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where we also used conservation of norm of u. Hence, if we choose Eq < S we can define 
h{a{t)), t G M, see Proposition 12. 1[ We then obtain (fTT|) where 

r{t) = u{t) - a{t)tPo - h{a{t)), (V^o, r{t)) = 0. 

The solution is now described by the scalar a{t) e C and r{t) G C{M.,H^). Their equations 
are obtained by projecting onto tpQ and its orthogonal complement in : 

= E{\a{t)\)a{t)+^{^o.2\^E\^r + ^lf + 2ijE\r\'' + ^Er^ + \r\h) (12) 



dt 

dv — 
i— = Hr + 7PJ2|?/^b| V + ?/;|f + 2?/;^^|rp + i^Er^ + \r\^r] 
at 

-lDh\ai^t) 2\4jE\^r + tPlf + 2tPE\r\^ + ^Pet^ + |r| V) (13) 



where we used the identities ipE = aipo + h{a) and _D/;,|a[a] = E{\a\)h{a) for a G C, \a\ < 6. 

In order to obtain the estimates for r{t), we analyze equation p3|) . In the next section we 
study its linear part: 

= Hz + ^P,[2\iPe?z + - 7D/i|,(j)(^o, 2\iPe?z + ^Iz)] 




Let us denote by f2(t, s)w the operator which associates to the function v the solution of 
the above equation: 

fi(t, s)v = z (14) 

The estimates that we need for this linear propagator are proved in the next section. 
Now, using Duhamel's principle (fT^ becomes 

r(t) =n(t,0)r(0)+ / n{t,s){-iPc[2^EHs)\^ + i)Er''{s) + \r\^r{s)] 
Jo 

—fDhlit) [(V^o, 2^E\r\^ + ^eT^ + \r\^r)]}ds (15) 

It is here where we differ essentially from the approach for the 1-d case |21 13 El and 3-d 
case j?! 123 I2H1 El- The right hand side of our equation contains only nonlinear terms in 
r. Hence, if we make the ansatz r{t) ~ (1 + t)~^^^ then the quadratic and cubic terms in 
(fTH|l decay like (1 + s)-^/^ respectively (1 + s)"^/^ Both are integrable functions in time, 
hence, via convolution estimates, the integral term on the right hand side decays like Q{t, 0). 
We have a chance of "closing" the ansatz provided f2(t,0) ~ (1 + 1)~^/^. Contrast this with 
the case in the above cited papers where a linear term in r is present on the right hand 
side. Same argument leads to a loss of 1/4 power decay in the linear term and requires 
Q(t, 0) = Q-'tHt ^ ^-^ _|„^^-i-'5^ S > 0. for closing. This turns out to be impossible in norms 
in 2-d, see (j24p . while the use of weighted norms, see (I23j) . for delocalized terms as the 
cubic term in psp. would require compensation via virial inequalities, see [21, which needs a 
much higher power nonlinearity than cubic in the delocalized terms^. 
The following Lemma makes the above heuristic argument rigorous: 

^Heuristically we arrived at quintic power nonlinearity, hence this technique may be apphcable to the 
quintic Schrodinger in 2-D but definitely not to the cubic one. 
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Lemma 3.1 There exists si > and 62 > such that if \\ < x >" il)E\\m < O'nd the initial 
condition r(0) satisfies 

max{||r(0)|U2,||r(0)||Hi} <£2 

the initial value problem / f73)) is globally well-posed in C(]R, L^^ fl fl L^), 6 < p < 00 and 
for allteR: 

l|r(t)|U. < 



(1 + \t\y-^/p 

C'2eolog(2+|t|) 



\r(t)\\L^ < 0360 



with Ci, C2, C3 independent of Eq and C2 depending on p > Q. 

Note that the Lemma finishes the proof of the theorem. Indeed, we now have two solutions 
of ((131), in C(M,H^) from classical well posedness theory and one in C(M, fl fl 
L^), P > 6 from the Lemma. Using uniqueness and the continuous embedding of in 
L^^ n n L^, we infer that the two solutions must coincide. Therefore, the time decaying 
estimates in the Lemma hold also for the solution. The L^, 2 < p < Q estimates in the 
theorem follow from interpolation: 



C'3C'2£olog(2+|t|) 



{i + \t\y-^/p 

It remains to prove the Lemma: 
Proof of Lemma 13.11 

Fix p > 6. We will show that (fTH|l has a solution by applying the contraction principle in the 
functional space 

Y = {f:R-.L\nL^nL'\ sup (1 + ||/(t) 1^. < 00 

t>o 

SUP \ WfmL. < 00,SUp||/(t)|U. < 00} 

t>o log(2+|t|) t>o 



endowed with the norm 



y = max{sup (1 + |t|)^'^||/(t)|U.^,sup \ ^/ ^ ||/(t)|U., sup ||/(t)|U.; 
t>o t>o log(2 + \t\) t>o 

To this extent we consider the operator defined on functions in Y as 

{Nu){t) = n{t,s)v+ [ n{t,T)-f{Pc[2'^E\u\'^ + 'ipEu'^ + \u\'^u] 
Jo 

-Dh\a(r){'ipO,'^'ipE\u\'^ + "ipEU^ + \u\\)}dT 

We will need some properties of the operator which are summarized in 
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Lemma 3.2 We have 

(i) The range of N is Y i.e. N : Y ^ Y is well defined, 
(a) There exists C > such that 

\\Nui - Nu2\\y < C{\\ui\\y + \\u2\\y + Wuify + ||M2||y)|ki - U2\\y 

In particular N is locally Lipschitz. 

Moreover, C = C{C,Cp,Cp^g'^) where the constants C,Cp,Cp g'^ are those from the linear 
estimates for Q{t, s) (see Theorems \4 1[ ^^in the next section). 

Proof of Lemma 13.21 Let us observe that it will suffice to show part (ii) and then using the 
fact that A^(0) = we will have part (i). Indeed, part (ii) will give us that for ui,U2 & Y we 
have Nui — Nu2 € Y . Taking ^2 = and since iV(0) = this will imply that Nui G Y . 
Thus, take ui,U2 € Y and consider the difference Nui — Nu2, which is 

{Nui - Nu2){t) = [ Q{t, r)-f{P,[24jE{\ui\ - \u2\){\ui\ + \u2\) + 4je{ui - U2){Ui + U2) 
Jo 

+ {Ui - M2)|MiP + (|Mi| - \u2\)iu2\Ui\ + M2IM2I)] 

-Dh\a{r){i^o,'2i^E{\ui\ - |n2|)(|ni| + |n2|) + - U2){ui + U2) + 

(ui - n2)|nip + (|ni| - |n2|)(u2|ni| + U2\u2\))}dT 

The L^^ estimate We can work under the less restrictive hypothesis: A < p < oo. Let If' 
be the dual of LP, i.e. 1/p' + 1/p = 1. We have 

\\Nu,-Nu2\\l2^< [ Mt,T)\\L2^L2^X 

Jo 

x|| 2iIje < X (|mi| - |m2|)(|mi| + \u2\) + ipE < X {ui -M2)(mi +M2) U^dr 

" V ' 

A 

Jo 

x|| {ui - n2)|ni|^ + (|ni| - |n2|)(n2|ni| + U2|'«2|) Wlp'^t 

" V ' ^ V ' 

Bi B2 

Jo 

x{\{'^jJo,2^|JE{\ul\ - |n2|)(|ni| + \u2\) + ^Pe{ui - U2){ui + U2) > | + 
^ V ' 

F 

+ \ <1pO, {Ui - U2)\ui\'^ + {\ui\ - \U2\){U2\UI\ + ^2 1 ^2 | ) ) | jc^T (16) 

" V ' 

G 

To estimate the term A we observe that 
II <a; >"^b(|mi| - |m2|)(|mi| + |m2|)|U2 < || < a: 7/;e|U.||mi - M2|Up|||mi| + |m2|||lp (17) 
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with i + ^ = i Then 

a p 2 



f \m,r)\\L._,L. A{r)dr 
Jo 



< I TT-T. . iM_ 2/0 , I. —■m'E<X>^\ui-U2\{\ui\ + \u2\)\\L^dT 



c_ 

(l + |t-r|)log^(2+|t-r|) 

<3^^ ;* log^(2+|r|) lllnil-ln^llly \\W\ + \u, 



(1 + |t - r|) log^(2 + |t - r|) (i + (i + 

\\Ul - U2\\y 



<3CCiC2(||ni||y+||n2|h 



(l + |t|)log^(2+|t|) 



where for the first inequahty we used Theorem I4.11 part (i). The constants are given by 
Ci = supi>o II < x ^i^lli. and C2 = supi>o(l + l^l)log^(2 + |t|) /J 



log"(2+|r|) 



dr 



(l + |t-rl)log^2+|t-r|) 



(l + |r|)'-p 



T < 00, because p > A. 



To estimate the cubic terms Bi, B2 we can not use the term ipE as before, and this is 
what forces us to work in the space. We have: 

\\{ui - n2)|uip||^p' < ||ni - n2||Lp||^ii||ic 

respectively 

\\{ui - n2)(n2|ni| + n2|n2|)||ip' < ||ni - n2||Lp||^i2||L"(||^ii||L- + lk2||L") 
with ^ + i = J^. Since 4 < p we have 1 < a < p. Therefore we can again interpolate: 

ll'^ijllL" ^ ||^i||j;^2 II II LP ;^ — 1)2, 

where - = + -. Combining these relations we obtain for : 

||(ni - n2)|nip||^p' < ||ni - n2||Lp||^ii||^2"^~^^||^ii||ip (18) 

respectively, for B2 ■ 

\\{U1 - U2)(^^2|^il| + ^i2|^^2|)||Lp' < ll^^l - '«2|Up||^^2||i2^||^i2||Lp(||^il|li2^||^il||LP + || ^^2 || ^2 || ii2 1| LP ) 

(19) 

with 

2(1-6) 2b 1 1 

— H h - = — . 

2 p p p' 

A consequence of this relation and of p < oo is: 

(l-^)(l + 26) = l + 2/p> 1 (20) 
which will play an essential role in what follows. 
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Thus, the estimate for the term containing Bi + B2 is 



||fi(t,r)Pc|Lp'^L2 \\Bi + B2\\lp' < Cp(||ni||^ + ||u2||y)||ni-n2||y X 

— fj 

^ log(2+|r|)(i+^^) 1 



where for the first inequahty we used Theorem 14.11 part (ii), inequahties ()18p. p9|) and the 
definition of the norm in Y. For the last inequahty we used the fact that (1 — -)(1 + 2b) > 1 

(see m) with Cs = sup,>o(l + \t\Y-^/^ £ '"^^^^^'-fr" , < ~- ' 

For estimating the term containing F we have 



l-^l < 3||V^o||L°°||V^i?||L-||^il - II LP (II '"ill LP + |k2||Lp) 



- H — = 1. Then, the term containing F is estimated as the term containing A with Ci 

replaced by C4 = sup^^p ||^/i|a(r)||L2 IIV'oIIl-IIV^eIIl"- 
We estimate G as 

IGI < 11^0 < a; ||lc||mi - MslL^^dkillip + ik2||ip + ||mi||lp||m2||lp) 

with - + i + - = 1. Then the term containing G is estimated as 

riior.^^ii ir^M^^o r Cll^o < ^ lU^ log'(2 + |r|) 

X (1 + It - r|)log^(2 + It -r|) -(1 + 1x1)3-/.^^ 

< sccsdi^iii?- + 11^211^) ., , l"!;,'?"^ — 



with ^5 = IIV^o < X SUp,>o(l + N) log'(2 + |t|) /o (i+|t-H)lo.^72|?J:|'!(l+k|)^-^/^ ^^ < ^ 

because p > 3. 



(l + |t|)log^(2+|t|) 

log'(2+|r|) 

(l + |t-T|) log2(2+|t-r|){l + |T|)3-6/P ' 
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The estimate: With p', q'q, q' given by Theorem 14. 2[ we have 

\\Nui- Nu2\\lp < [ Mt,T)\\L2^LP X 

Jo 

x|| 2ijjE < X {\ui\ - |n2|)(|ni| + |U2|) + i>E < x {ui - U2){ui + U2) lUadr 



A 



+ f\m.r) 

Jo 



'LP'nL'''onLi'-*LP ^ 



x|| (ui - n2)|ni| + (|ui| - |n2|)(n2| 

Bi B2 

+ [ \\Q{t,r)\\Ll-.Lp\\Dh\a(r)\\Ll^ 

Jo 

x{\{^Q,2tpE{\Ul\ - |^i2|)(|^^l| + |^i2|) + i'E{Ul - U2){Ui + ^2))! + 

V ' 

F 

+ Ktpo, {Ui - U2)\Ui\'^ + (|mi| - \U2\){U2\UI\ + M2 1^2 | ) ) | (21) 

^ V ' 

G 

The term A can be treated exactly as before and for the ^2 term we use Theorem \4. II part 
(Hi). Since 1 < p', q'^, q' < 4/3, we can estimate the -81,-82 terms in each of the norms 
LP , L'^°, L'^ as we did above for their norm only. For Q we use Theorem 14 ■2[ part (Hi). 
The terms F and G are also treated as in the previous case. The convolution integrals in 
(12 ip will all decay like (1 + except the second one which will have a logarithmic 

correction dominated by log(2 + \t\). 

The estimate: We have 

\\Nui-Nu2\\l^ < f Mt,T)\\L2^L2X 

Jo 

X II 2iIje <x>" {\ui\ - |n2|)(|ni| + |n2|) + ^/'s < a; [ui - U2){ui + U2) Wi'^dr 

V ' 

A 



X|| (mi - M2)|mi|^+ (|mi| - |m2|)(m2|mi| + M2 |m2 | )J| Lp'nL2lir 



Bi B2 



+ [ \\n{t,T)\\L2^L2\\Dh\a(r)\\Ll^ 

Jo 

x{|(^o,2^i=;(|Mi| - |m2|)(|mi| + IM2I) +fpE{ui - M2)(Mi + U2))\ + 

" V ' 

F 

+ \{ipo, {ui - n2)|ni|^ + (|ni| - |n2|)(n2|ni| + U2\u2\))\}dT (22) 

V ' 

G 

We estimate the term A as in ()17p while the estimates in L^' for Bi, B2 are as in p8|) and 
p9|). For their estimate in L?' norm we use 

||(mi - U2)\Ui\^\\l^ < \\ui - n2||Lp||^/i||i<:< 
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respectively 

\\{ui - U2){u2\Ui\ + U2\u2\)\\l2 < \\ui - ^2 || LP 1^2 ( H^l || + ||m2||l") 

with ^ + i = i. Since 6 < p we have A < a < p. Therefore we can again interpolate: 



where - = Kr + -• CombininEf these relations we obtain for i?i : 

11^ M |2|| ^ II II II l|2{l-fe)|| ||26 

\\{ui-U2)\ui\ ||l2 < - M2||lp||mi||^2 ||mi||lp 

respectively, for B2 : 

||(M1 - U2){u2\Ui\ + U2\u2\)\\l2 < \\ui - ^2 Hlp || ^2 || ^2 1| ^2 || ( || ^1 1| ^2 1| Ml || + || ^2 || ^2 || ^2 1| ; 

with 

2(1 -b) 2b 1 1 

— H h - = -. 

2 p p 2 

A consequence of this relation is: 

(l--)(l + 26) = 2 
p 

Using now the definition of the norm in Y we will have: 

11^1 + B2\\l < \\Ui -M2||y(||Ml||y + ||m||2 



2.1,., I, n,., 1,2 , ||.,||2^1og'^/^"-'^(2+|t|) 



{i + \t\y 

The previous estimates for F and G suffice here as well. 

Recalling from Theorem \4.1\ part (in) and Theorem \4.2\ part (ii), that r)||i2^/^2 
and t)||^p'|-|2.2^l2 are bounded, and combining with the estimates above, as well as 

taking into account the definition of the functional space Y we have that 

llA^Mi - Nu2\\l^ < C\\ui - M2||y[C'6(||Mi||y + ||M2||y) + CriWuify + ||M2||y)] 

with Ce = supj>o /J j^^^^dr < oo and C7 = supt>o /g (YTmF™^^^^ ^ ^■ 
This finishes the proof of Lemma \I?A □ 
We can continue now with the proof of Lemma ^^ Let 

b = minjl — e, — ^| 
AC 

where e > arbitrary and C is given by Lemma IH.2I Consider the ball of radius b centered 
at 

B = B{0,b) 

Using Lemma^^ part (ii) with U2 = we have, for ui E B 

\\Nui\\y < C'dlwilly + ||ni||y)||ni||y <C2b-b<b 
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which means that the ball B is invariant under the action of the operator A^. 
Also, using again Lemma ^TA part (ii), we have, for mi,M2 G B 

\\Nui - Nu2\\y < C{\\ui\W + ||m2||y + \\ui\\y + ||M2||y)||'"l - U2\\y 

< CAh\\ui - U2\\y < (1 - e)||ni - n2||y 

which shows that : i3 ^ i3 is a contraction. This finishes the proof of Lemma 13.11 and of 
Theorem \'d.l\ □ 



4 Linear Estimates 

Consider the linear Schrodinger equation with a potential in two space dimensions: 

zf = {-A + V{x))u 

m(0) = Mo. 

It is known that if V satisfies hypothesis (Hl)(i) and (ii) then the radiative part of the solution, 
i.e. its projection onto the continuous spectrum oi H = —A + V, satisfies the estimates: 

lle-'«P..„|U,. < gM (i^|,|).4.p^|,|) ll"olUl (23) 
for some constant Cm > independent of Uq and t G M, see ^3 Theorem 7.6], and 

||e-^^*PeWo||L. < ^^IkolliP' (24) 

for some constant Cp > depending only on p > 2 and p' given by p'~^ + = 1. The case 
p = cx) in (|24p is proven in [221 • The conservation of the norm, see Corollary 4.3.3], 
gives the p = 2 case: 

\\e"'^^PcUo\\L2 = \\uo\\l2. 



The general result ()24|) follows from Riesz-Thorin interpolation. 

We would like to extend this estimates to the linearized dynamics around the center 
manifold. In other words we consider the linear equation, with initial data at time s, 

'z| = (-A + V{x))z + ^P42\tljE{t)\^z + ijUt)z + Dh\ait){{^Po, 2\iljE\^Z + ^1^))] 
z{s) = V. 

Note that this is a nonautonomous problem as the bound state ipE around which we linearize 
may change with time. 

By Duhamel's principle we have: 



z{t) = e-'^^'-'^ PXs) - i j e-'"^'-^hP^[2\^E\^z + ^Iz + 

Dh\a(r){{M'^E?z + ^lz))]dr (25) 
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As in fll4|) we denote 

Q{t,s)v=^z{t). (26) 

In the next two theorems we will extend estimates of type (|23 p - (|24p to the operator Q{t, s) 
relying on the fact that ipE(t) is small. It would be useful to find sufficient conditions under 
which our results generalize to large bound states. Such conditions have been obtained in one 
or three space dimensions, see [21 ESI 123 EI , unfortunately their techniques cannot be applied 
in the two space dimension case. 

We start with estimates in weighted spaces: 

Theorem 4.1 There exists ei > such that if \\ < x >^ '^PeWh'^ < then there exist 
constants C, Cp > with the property that for any t, s G M the following hold: 

C 

(i) \\VLit, s)\\t2^t2 < ; TT, ; 

(J 

(a) \\^{t, s)\\tp' < — 2") for ^-^y oo> p>2 where = 1 

— a \ I I 1 — — 

\t — s\ p 

C 

(iii) \\^{t,s)\\Li^LP < for any p > 2 

|t — s\ p 

Before proving the theorem let us remark that (i) is a. generalization of while (ii) and 
(Hi) are a mixture between (j^Hj) and (j23). We have used all these estimates in the previous 
section. They are consequences of contraction principles applied to (j23) and involve estimates 
for convolution operators based on ()23p and (j211)- It will prove much more difficult to remove 
the weights from the estimates (ii) and (Hi), see Theorem 14.21 

Proof of Theorem 14.11 

Fix s e M. 

(i) By definition (see (j^ ). we have Q{t, s)v = z{t) where z{t) satisfies equation (j^ . We 
are going to prove the estimate by showing that the nonlinear equation can be solved 
via contraction principle argument in an appropriate functional space. To this extent let us 
consider the functional space 

Xi := {z G C(M, L2^(M'))| sup(l + \t - s\) log2(2 + \t - s\)\\z{t)\\L^ < oo} 

tm 

endowed with the norm 

\\z\\x, := sup{(l + \t- s\) \og\2 +\t- s\)\\z{t)\\L2 } < oo 

Note that the inhomogeneous term in (j2Hl): 
satisfies zq G Xi and 

\\zo\\x, < CmMli (27) 
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because of 

We collect the z dependent part of the right hand side of in a linear operator : 
Xi — > Xi, 

J s 

In what follows we will show that L is a well defined bounded operator from Xi to Xi 
whose operator norm can be made less or equal to 1/2 by choosing Si in the hypothesis 
sufficiently small. Consequently Id — L is invertible and the solution of the equation ()25|) can 
be written as z = {Id — L)~^zo. In particular 

\\z\\x, < (1- \\L\\y^\\zo\\x, <2||2;o|Ui 

which, in combination with the definition of Q, the definition of the norm in Xi and estimate 
(|27j), finishes the proof of (i). 

It remains to prove that L is a well defined bounded operator from Xi to Xi whose 
operator norm can be made less than 1/2 by choosing ei in the hypothesis sufficiently small. 
We have the following estimates: 

\\L{s)z{t)\y^ < 1^ ||e-^(*-^)Pe|L._.,.^ ■ [3|||^^p(r).(r)|U. 

< r We^'^'^'-^^PcU^Ll^ ■ [3mE\\T)z{T)hl 
J s 

+ ||L'/l|a(r)||c^L2||^o|U2 3||^|||L-||2(r)||i2_J 

On the other hand 

III^Er^^lUa < \\4liJ < X >^'' ll^i^riU-, and II < X ^ij||ioo < el (28) 

where the last inequality holds because of the Sobolev imbedding if^(M^) C L°°(M^) and of 
the inequality 

II < X ^e\\h^ < £i- 

Also 

||-D/i|a(r) II < C^l, as k(''")| < 

Using the last three relations, as well as the estimate (j23p and the fact that z E Xi we obtain 
that 

\\L{s)\\x,^x, <£iSup[(l + |t-s|)log2(2+|t-s|) X 

t>0 



1 

1 + |t - r| log'(2 + |t - r|) (1 + |r - s|) log'(2 + |r - s|) ^ 



X / ,^ ^ ,^ ^ ■ 7^ ^ , -dr < 2Ci£i (29) 
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Indeed, in order to prove the above we will split X into A + B where 



. t + s 



~ 1 1 

A = I 7i 7i dr 

'l + \t- r|) log'(2 + \t- r|) (1 + |r - s\) \og\2 + |r - s|) 



for which we have the bound 

\A\ < 



t + a 

1 . r~ 1 



1 + I VI) log'(2 + I VI) Js (1 + k - s\) \og\2 + \t-s\) 

1 



;i + IVI)iog^(2 + lVI) 

Observing that A = B and using the last estimate in ()29|) we obtain that 

Ili^lUi^Xi < Cisi < 1/2 

for £1 small enough. 

(a) By the definition of Q it is sufficient to prove that the solution of (j25p satisfies 

C 

\\z{t)\\L2_^ < ^^tII^IIlp'' for all oo >p > 2 where p'^^ + = 1. (30) 

|t — s| p 

We will use a similar functional analytic argument as in the proof of (i). Fix p,2 < p < oo 
and assume f G , p'~^ + p~^ = 1. We will work in the following functional space: 

X2 := {z G C(M,L^^(M2)|sup||^(t)||i2 |t - < 00} 

endowed with the norm 

1— - 

lkll^2 •= sup ||z(t)||j;^2 |t — s| P < 00. 

Using the fact that L^^. continuously and the estimate ()24|) we have e'"^^^^"^^ PcV G 

X2. In addition, for L defined in the proof of (i), we have 



1 

sup|t — s| ^p\\L{s)z{t)\\i^2_^ < 



t>0 



sup It - s^-p f ||e-^(*--)Pe|U|^L^_^ • mE?{r)z{r) + 
+V'|(r)z(r)|U2 + ||P,D/i|,(,)||c_L2KV^o,2|V^i5|'z(s) + ^/^|z(s))|rfr 



< sup It - .1-^ / ' h-^-^h-^h;^-— hLo. ^ ^31) 



:i + |t - r|)|r - s\'^p \og\2 + \t - r|) 



Using now the bounds (j28|l in (j3T| . for ei small enough, we obtain that the norm of L{s) is 
less or equal to 1/2, i.e. the operator Id — L{t, s) is invertible, which, as in the proof of (i), 
finishes the proof of estimate (ii). 
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(in) We already know from part (i) that equation ()25p has a unique solution in L^^. 
provided v E L"^. We are going to show that the right hand side of is in L^. Indeed 

\\e-^^i^-^^PMs)\\Lr < , \ . Hs)h, < \ A v{s)\\,. (32) 

\t — s\ p \t — s\ p 

where the Cp's in the two inequalities are different, for the first inequality we used fl24|) while 
for the second we used the continuous embedding , 1 < < 2. For the remaining 

terms we combine ()24|) with < oo obtained in part (i): 



* 2a 



J s 



J s \t — t\ p 

< / dr < ^ (33) 



s\ p 



with i + 1 = ^. 

Similarly, we have 

< f ^^^\{^O.mE?z{s)+,ljl-z{s))\dT 

J s \t — t\ p 

pt (J 

< / ^jzi\\\^o\\l'A\ < X >" i^l\\L°-\\ < X z\\l2 

Js \t — t\ p 

/* 1 1 c 
2 o dr < (34) 
_ |t_^|i-f (l + |r-.|)(log2(2 + |r-s|) -\t-s\'-p 

Plugging (jS21)-(|nil) into ^ we get: 

C 

\\z{t)\\LP < ^tziIMlI 

\t — s\ p 

which by the definition Q{t, s) = z{t) finishes the proof of part (in). □ 

The next step is to obtain estimates for Q{t, s) in unweighted spaces. They are needed 
for controlling the cubic term in the operator of the previous section. 

Theorem 4.2 Assume that || < x "ipEWn^ < (where ei is the one used in Theorem \4.1\ ). 
Then for all t, s eM. the following estimates hold: 

.Ml ^ Cp^q>\og{2 + \t-s\) 

{l + \t-s\f-p 
for all p, q', 2 < p < oo, 1 < g' < 2, + p"^ = 1; 



17 



4 

(ii)\\n{t,s)\\^,^^,'^^^, < C,., for all q'o, 1< q'^ < 
(in) for fixed po > and 1 < gg < 4/3 and for any 2 < p < pq 

l-2/p 

C„„' log(2 + \t-s\)^=^ 



\t — s\ P 



where 



with 



i-2/p 1 e i-e 

6 = f-, i.e. - = \ . 

1 - 2/po P Po 2 

Note that (Hi) is similar to the standard estimate for Schrodinger operators (|24jl except 
for the logarithmic correction and a smaller domain of definition. We will obtain it by inter- 
polation from (i) and (ii). The proof oi (i) will rely on a fixed point technique for equation 
(jHUI) while the proof of (ii) will rely on Strichartz inequalities. 

It turns out that we need to regularize (j^^ in order to obtain (i) and (ii). The inhomo- 
geneous term has a nonintegrable singularity aX t = s when estimated in : 



-iH{t-s) p^n^^ < \t-s\-^\ 



Using estimates with integrable singularities at t = s, for example in L^, p < oo see 
would lead to a slower time decay m (i) and eventually will make it impossible to close the 
estimates for the operator in the previous section. We avoid this by defining: 

W{t) =^ [z{t) - e-'^(*-^)PXs) (35) 
which, by plugging in (^3]) . will satisfy the following "regularized" equation: 



W{t) = -I 



I f e~''"^'~^hPc[2\Mr)\'e-'''^^-'^PMs) + 4jUr)e'''^^~'^ Pcv{s)]dT 

J s 



~m 



J s 



— I 



fit) 

t 



V ' 

m 
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Some other new notations are necessary for the sake of easy reference. We will denote by 
T{t,s) the operator which associates to the initial data at time s, v, the function W{t), so 
that 

T{t, s)v =^ W{t) (37) 
which will be related to the operator Q{t, s) = z{t) (see fll4j) ) by 

n{t, s) = T{t, s) + e-'"^'-'^Pc. (38) 

For T we can not only extend the estimate in Theorem 14.11 (ii) to the case p = oo but also 
obtain a nonsingular version of it: 



Lemma 4.1 Assume that || < x i^EWn^ < (where e\ is the one used in Theorem \4-.l\} . 
Then for each 1 < q' <2 there exists the constant Cq> > 0, Cq/ ^ oo as q' ^ 1, such that for 
all t, s eM. we have: 



Cq' 



t 



Proof of the Lemma: Fix g', 1 < g' < 2. Consider equation ()3fi|l with s G M arbitrary 
and V G L^ClL'^ . We are going to show that ()36|) has a unique solution in C(M, L^^^) satisfying: 

^ l^^''_^| "^ax{||t;|Ui, \\v\\^,,} 

which will be equivalent to the conclusion of the Lemma via the definition of T ()37p . 

Let us observe that it suffices to prove this estimate only for the forcing term f[t) + f{t) 
because then we will be able to do the contraction principle in the functional space (in time 
and space) in which /(t) + /(t) will be, and thus obtain the same decay for W as for f(t) + f{t). 

Indeed, this time we will consider the functional space 

X3 := {u G C(R,L^^(R2)|sup||M(t)||i2 (1 + |t - s|) < 00} 

tm 

endowed with the norm 

||m||x3 := sup{||M(t)||i2_ (1 + |t - s\)} < 00 

We have 

sup(l + |t-s|)||L(s)n(t)|U2 < 
t>o 

sup(l + \t- s\) [ ||e-'^(*-")p,|L2_i2 X [2|| <x>'' < x < x n(r)|U2 

+ ||L'/i||i2||V'o|U2|l <x>'' ^plh^lHr - s)\\LlJdT 

< sup(l + \t~s\) —— "L^i r, < Cssi (39) 

t>o Js (1 + |t - r|)(l + |r - s|) log (2 + |r - s|) 

Using the bounds (j28|l in (j39|) we obtain that for ei small enough the norm of L{t, s) in 
X is less then one, i.e. the operator Id — L{t, s) is invertible. 
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We need now to estimate f(t) + f(t): 



(where we used the estimate 

We will split now (pOj) into two parts to be estimated differently: 



cs+l pt 

II/ + /IIl^„< / ...+ / ••• (41) 

Js+1 



1 XX 

Then, we have: 



' '"A (l + |t-r|)log2(2 + |t-r|) 

< ^ ^-4^- . ^ / \\e-'''^^-'^Prv(s)\\L. ■ III^Ep < X lU. < 

V ' 

< fixed constant 

< —0 \\vis)\\rA^—Y"^dT 

- il + \t-s-l\)log\2+\t-s-l\)Js " ^ V-s^ 



il + \t-s-l\)\og'i2 + \t-s-l\)J, ^ 



/ C'io|I^(^)IIl.' 1 IL.^„MI 

- (l + |t-.-l|)log^(2+|t-.-l|) - + 



with i + i = i and i + 4 = 1. 
For the second integral we have: 

Js+1 (l + |t-r|)log2(2 + |t-r|) 

- (H-|t-r|)log^(2+|t-.|) ■ 

Let us observe that the last two estimates are for the case when t > s + 1. Ifs<t<s + 1 
we have 

"^ + ^"^--1 (l + |t-r|)log^(2+|t-r|) 

- (l + |t-r|)(log2(2+|t-r|)) 

<Ci3 r(^)'-'cir||t;(.)||,,<C||t;(.)|L,, 

with i + 1 = 1 

a q 2 
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Combining the last three estimates we get the lemma. □ 
We can now proceed with the proof of Theorem 14.21 
Proof of Theorem 14.21 

(i) Because of estimate (j^^ and relation (jHHj) it suffices to prove (i) for T{t, s). 

Consider equation (j36|) with arbitrary s G M and v & (1 L"^ . In the previous Lemma we 
showed that the solution W{t) G L'^„- Now we show that it is actually in for all 2 < j» < oo. 
Fix such a p. Then: 



\\wmL,<\\f{t) + m\\Lr+ I iie"^^(*-^)p,ii^,._i,(2+iiD/iU(.)iUpiiV'oiiL.) 



<\\fit)+m\\ 



LP 



2C 



{t-ry-v 



< X lU^II < a; W{T)\y (42) 



The estimate for f{t) + f{t) is similar , but this time the term 
controlled for s + 1 < r < t by 



eWlp 



I e 



C 
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\t-s\ 



|2g-i/f(r-s)p^^||^^, is 



and for s<r<s + lby 



C 
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r — s 1 



where a^^ + q^^ = p'~^ and + q'~^ = 1. 

Using now the previous Lemma to estimate the term || < x >~'^ W{t)\\l2 and replacing 
in (jl21) we get: 



I|w^WIIl.< 



Ci7log(l+|t-s|) 



2 ■max{||i;||ii, ||?;||^,/} 
[l + \t-s\Y-p 



with 1 < q' < 2 which is equivalent to 



Ci7l0g(l + |t-s| 

(l + |t-s|)'~i 



for all 2 < p < oo and 1 < q' < 2. This finishes the proof for part (i). 
(a) Recalling the equation for W flHUj) . let us observe that we have 



(43) 



e-^^(*--)p,(|^^|2iy(r) + i>lW{T))dr\\L2 < Cs{ / Ue\'W{ 



<C, 



<x>'' f , , WWir) < X \\\dr)T < Cssl 



2y 



7'(1--) 



dr 



{l + \r-s\) 

<CMk (44) 
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where for the first inequahty we used the Strichartz estimate 

(T/)(t) = f e~'^^'~^^f{T)dT : L"''{0,T-LP') ^ L°°{0,T-L^) 

J s 

with (7, p) with 7 > 2 . For the second inequahty we used Holder's inequality and for the 
third one we used ()30|) combined with (jH^j) and (j211). Finally the last inequality holds when 
7'(1 — ^) > 1 which happens for qo > 2j > 4. 



2_ 

QO 

Also, we have the estimates 



<Cs{f p/iU(.)(V^o,2|^i^riy(r) + ^|iy(r))||^>r)7 

J s 

<e,C,,if\{^o,2\^E\'W{T)+^lW{TW'y 

J s 

(1 + |r-s|)^'"i^' 



< ^iC2i( / :jrT^dTy\\v(s)\\^,, < e^CM^,, (45) 



where for the first inequality we used Strichartz estimate as before and for the second inequal- 
ity we use the fact that Dh\a{T) is bounded in and thus in any and its norm is small. 
For the fourth inequality we used the fact that ||'?/'o||l2 and |||'?/'_eP < x H^oo are bounded 
and small. Finally the last inequality holds, as before, for go > 27 > 4. 

For f{t) + f{t) we'll need to estimate differently the short time behavior and the long time 
behavior, namely: 

+1 rt 



...+ / 

J s- 



+1 

X II 



We have: 

s+l 



\TW < 2|| e-^^(*-^)7Pe[|^i5|'e-^^(^-^)p,t;(s) + D/i((^o, 2\ijE\'e-'" 

J s 
J s 

.s+l ^ 

^ ^-25 / ;izx^^II^IIl^ ^ C2e\\v\\^,> 

J s [T — S) ''0 

(where we used the fact that the operator e~^^^ preserves the L? norm, and ^ + ^ = \)- One 
estimates similarly the terms containing ?/^|^ instead of IV'eP- 
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We continue by estimating XX: 



Js+l 

Js+l 

+( r ||V^olUHIlV'i.riUH|e-^''^^"'^^'||2'.orfr)i^ 



< C27( /* ||e-^^(^-^)Pc^'||^:)^ < C,,{ f l^i^^dT)^\v\\,, < C,4vu 



s+i Js+l [r — s) "JO 

where for the first inequahty we used the fact that the L/^ norm is preserved by the operator 
e~^^^Pc- For the second inequahty we used the Strichartz estimate 

(Tf) it) = I e-'"'^'-^^f{T)dT : L'^o (0, T; L"') ^ L°°(0, T; L^) 

J s 

for the f{t) term. For the f{t) term we used similarly the same Strichartz estimate, the fact 
that llD/ijl^g/j is bounded (as it is in any norm), and we estimated the scalar product by 

the product ||'?/'o||L2|llV^-Bp||L"||e~*''^*-'^~'^-*f (s)||i9o. For the third inequality we used Holder's 
inequality and the fact that the ||'?/;£;p||^^i, || iV'sPlU"; ll'^olU^ < C,Vt ( where ^ ~ qo ~^ ^ 
1 = ^ + Finally the last inequality holds because — > 1, as go > f • 

Let us observe that we assumed that t > s + 1. Ifs<t<s + 1, only the estimate for X 
will suffice, where the upper limit of integration s + 1 should be replaced by t. 

Combining the estimates for X, XX, (j44p and (j45p we have that W{t) is uniformly bounded 
in which, by (|?7|). implies 

^)1L^^L2 < for all t, seR. (46) 

Using now and with p = p' = 2 we obtain (ii). □ 
(in) We start from ()43|): 

C„„, log(2 + \t- s\) 
(l + |t-s|) P 

and (gni): 

\\nt,s)\\^,,^^,<c,,, i<go<^- 

We can now use the Riesz-Thorin interpolation between the spaces fl L^o and L'^'o as 
starting spaces and between and as arrival spaces to get the claimed estimate. Indeed, 
it suffices to take as in the statement 9 = t—M^, and use it with the above two relations to 
get 

WTU M\ ^ C,,,^log(2 + |t-.|)^ 

W-^ V^^^JllLi'nLp'nL'^o^LP - u 

\t — s\ P 
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with 1 < < |,p > 2 and 




_ ^ l-2/p 

PO ' 2 l-2/po 



Using now ()H8|1 and we obtain the claimed estimate for Q{t, s). □ 



5 Conclusions. 

We have estabhshed that the solution starting from small and localized initial data will 
approach, as t ^ ±oo, the center manifold formed by the nonlinear bound states (solitary 
waves). However we have not been able to decide whether the solution will approach exactly 
one solitary wave as in the 1-d and 3-d case, see for example |3jnni- Here is the main reason: 
The long time dynamics on the center manifold is given by the equation p2|) . Since 

t 



a(±oo) — a(0) = lim [ —dt 
t~*±oo Jq dt 





the existence of an asymptotic limit at t = ±oo is equivalent to the integrability of the 
right hand side of (fT^ on [0, +00) respectively (— oo,0]. The terms containing and 
are absolutely integrable because they are dominated by (1 + \t\)'^^'^^P~^\ respectively (1 + 
|^|)3(2/p-i)^ which are integrable on M for p > 4. However, the linear terms in r do not decay 
fast enough to be absolutely integrable. It is possible though that a combination of decay 
and oscillatory cancellations would render it integrable. We think that is only a matter of 
time until a suitable treatment of this term is found. Note that, in the 1-d and 3-d cases, 
the linear terms in r were absolutely integrable in time, see for example [HE]. But these 
estimates relied on the integrable decay in time of the Schrodinger operator in L°° norm in 
3-d, respectively on the large power nonlinearity to compensate for the linear growth in time 
introduced by virial type estimates in 1-d. None would work for our cubic NLS in 2-d. 

The situation is even more complex and possible more interesting when the center manifold 
has more than one branch (more than one connected component). For simplicity, consider the 
case when hypothesis (HI) part (iii) is relaxed to allow for two, simple, negative eigenvalues 
Eq < El with corresponding normalized eigenvectors tpo, i^i- In this case the center manifold 
has two branches ipE^ = (ij'ipj + hj{aj), j = 0,1, each bifurcating from one eigenvector as 
described in Section |21 The decomposition into the evolution on the center manifold and the 
one away from it will now be: 

1 

u{t,x) = ^{aj{t)^j{x) + hj{aj{t)))+rm{t,x) 

j=0 

' V ' 

1pCM{t) 

The equation for r^lt) remains essentially the same as ()13|) in Sectional with ipE replaced by 
ipCM and the differential of h replaced by the sum of the differentials of hj, j = 0,1. However, 
one has to add to the right hand side of ()13|1 the projection onto the continuous spectrum of 
the interaction term between the branches: 

2| V'Eo I Vi^i + ^loV'iJi + SV'Eo iV'Ei r + i^Eo^k (47) 
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In principle one could use our techniques and obtain a decay in time for rm{t), hence collapse 
on the center manifold, provided one makes the ansatz that the term above, or at least 
its projection onto the continuous spectrum, decays in time. Such an ansatz needs to be 
supported by the analysis of the motion on the center manifold given now by a system of 
two ODE's, one for oq and one for ai. Each of the equations will be similar to (jl2p but the 
projection of (jT7|) onto ipO: respectively ipi^ has to be added to the right hand side. Note 
that, in the 3-d case, under the additional assumption 2Ei — Eq > 0, it has been shown that 
the evolution approaches asymptotically a ground state (a periodic solution on the branch 
bifurcating from ipo) except when the initial data is on a finite dimensional manifold near 
the excited state branch (the one bifurcating from ipi), see j2Zl EHl EOl EI] • But the authors' 
analysis relies heavily on the much better dispersive estimates for Schrodinger operators in 
3-d compared to 2-d. The 2-d case remains open. 

Returning now to the case of one branch center manifold in 2-d, an important question is 
whether its stability persists under time dependent perturbations. In jHj we showed that this 
is not the case in 3-d. The slower decay in time of the Schrodinger operator in 2-d compared 
to 3-d prevents us, yet again, from extending the technique in to the 2-d setting. 
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